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One-sided Levy stable distributions 

Jung Hun Han [j] 
Abstract 

In this paper, we show new representations of one-sided Levy stable distributions for 



irrational Levy indices of the type ( 1 1 1 which are not covered in [8j : for rational Levy 
indices. Furthermore, other equivalent representations for a distribution of a rational Levy 
index is described. We also give a simplest proof for the formulae which cover the cases for 
rational Levy indices. Finally we introduce the concepts of Levy smashing and Levy-smashed 
gamma stochastic processes. 

^ It was written in August, 10, 2010. 

^ 1 Introduction and preliminary results 

In |8] , Penson and Gorska obtained a general form of one-sided Levy stable distributions expressed 
as a Meijer G- function for rational Levy indices by putting v = and a = 1 in the formula 2.2.1.19 
in vol. 5 of [TO]. In [2J, the role of Mathai transformation in the theory of fractional calculus, 
which connects ordinary integral to fractional integral through their kernels, is described and a- 



level space is defined. Furthermore it is insisted that fractional integral and derivative preserve 

r(— — — ) 

the Levy structure defined in [2j. The Levy structure is nothing but the integrand ar ^_^ of 
t-h the if-function representation of the Levy distribution with a known as the Levy index which lies 

between and 1. Hence the case of simple irrational Levy indices of the type ( ? J 1 can be handled. 

£j In this paper, we provide formulae for one-sided Levy distribution of irrational Levy index for 

< a < 1 and some techniques to generate infinitely many new representations of one-sided Levy 
distribution. Furthermore, we introduce the concept of Levy smashing as a consequence of Levy 
effect on the family of gamma density functions and Levy-smashed stochastic processes. 

X 

We will use the following integral representation of the gamma function: 



poo 

F(z) = p z t z - l e~ pt dt, *K(p) > 0,9t(z) > 
Jo 

= lim — — — — — - — ■ — -, z ^ 0,1,2,3, ••• . 

n->oo z\Z + 1) ■ ■ ■ \Z + Tl) 

Pochammer symbol is defined as 

(6)* = 6(6 + 1)- .-(6 + A;-l), (6) = 1, 6^0 
T(a + k) , 

= {ajk = — ttt~^ — whenever the gammas exist. 
r(a) 
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1 



For if-function representations and their convergence conditions, consult with [6J, [?]: 



Trm,n 
PA 



I (ai,Ax),(o2,j42),— ,(a p ,A p ) 
l(6l,Bi),((62,B 2 ),-,((6 9 ,B g ) 



iUT=i £gj + B r s)HUU r ( x - a j - jig)) 
2« 71 {n?= m+1 r (i - 6, - ^)}{nt„ + i r (a, + 



z s ds. 



Generalized Wright function, which is well explained in [7] and which is a particular case of the 
^/-function, is 



I (ai,^4i),(o2,j42),— ,(a p ,A p ) 
l(6i,Bi),(6 2) B 2 ),-,(6 g ,B,) 



where a h b 3 G C and A h Bj G R: A ^ 0, ^ ^ 0, % = 1, • • • ,p,j = 1, • • • , g; £J =1 5, - £f =1 A, > 
— 1 for absolute convergence. 

Hypergeometric series, which is well explained in [5] and which is a particular case of the H- 
function, is 



,F q (at, ■■■ ,a p ;bx,-- - ,b q ;z) = 



n=0 



n?=i(^) n\' 



which is absolutely convergent for all z in C if p < q. 
Definition 1.1 Levy jump function is defined as follows 



J, 3 < n, 3,™ e Z" 

3+1 
9 : 



j > n, j,n G Z+. 



Lemma 1.1 [ Stirling asymptotic formula ] f^j/ 
For \z\ — > oo and a a bounded quantity, 



T{z + a) « (2vr) i /^ 



l/2 r ^+a-l/2 e -z 



Lemma 1.2 



iq-jp 



pq 



(_i)« r (j. 
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,) 
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.) 





pq 



and 



TTm,n 



where v,n,p,q,i,j G Z + and > 0. 



k H£f 



k I (o.p,kA p ) 
\{b q ,kB q ) 



Proof. Use the formula in [B], p44, + 1 — u) 
complete the proof. 



(-i)"r(/3+i) 



and the properties in [7j, which 
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2 Levy stable distribution of Levy index 



Let 



fa(x) 



^x- Sl d Sl , l>a> Re{s x ) > 0. 



2iri J L aT(l — si) 



(2.1) 



ar°i-s) a PP ears firstly in in the literature. Then is the well known Levy density function 

having the Laplace transform e~* Q . 



Theorem 2.1 Let < a < 1 and p,q,l\,l2 be arbitrary positive integers such that < a 



E q ) h < 1, p<q andl= (|V' 



-i 



/// is noi a positive integer and l\ ^ 1%, then 



IJpq 



IJpq 



^ iW),(a,0<->(^>0 



x(2tt) 



2_£ ^Zp^ I 9- 



-1 =*p 



xpy i(i,D,(|,i),-,(^,i) 



q-plq i\ ( (p-i)q—ylq 



x pl q q\l-l-l){i^-l)...^ 



+ 



IJpq 



x(2n) 



q-p 

2 



9-1 

E 



Ji R 

pr i 

x iq 3 



^5 ?l (^ ti -- , )(^,-l)-(fc^,-l) 



If I = 1 and l\ — h, then 
/(a,i)0) = 



_ 779- 1 . 

£ p— !)<?— 1 



pg 



x(2tt) 

X p-\F q _2 



x(2tt) 

q-1 

2 j=i \x q q 3 



3- 

pl q 



h. _ 


i 




g 


q 


h 






1 / «2<3- 


-jp 


) 


V P9 




(p-i)a-ip N 


pq ) 



If l belongs to the set {2, 3, 4, • • • } and l\ 7^ I2, then 



1 yqi-iJ V q iqiq-i 



f(<*,l)( x ) 



Ipql 



g-1,0 



III- 1 

fn \ g+l-P-t p+l-2,q-l 
X{27T) 2 



•^09 lf|,i),-,(^i,i) 



ar'g 



x(2tt 



pg7 ( P >lq ^ 



j V X J 9 gJ 



2=i- A 
9 LsJo-l' 



(P-I)g-J'pi 

pq ' 



(2.2) 



(2.3) 



(2.4) 
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Proof. Here we use a transformation 1 — si — as, then 



fa(x) 



T ^-x as - 1 ds, ->Re(s)> 



1 

a 



1 — a 



> 0. 



(2.5) 



2m J L T(as) a a 

Now, by using the Gauss-Legendre formula (Multiplicative formula) for gamma function, namely 

T(mz) = (2n) 1 -^m mz - 1 2T(z)T(z + -)---T(z+ r ^—^), m = l,2,--- , 

m m 



we have 
T(as) = T 



h 
p\'i 



Is 



ls\ „ ( Is 1 



v[p-\ = (2vr)^/T-5 r _ r - + - ---r - + 



q p 



Is p — 1 



Apply to the integrand in (2.5), then 



1 1 



r(s) x 



« 2™ y jr ( 27r ) VpPT-5r(^)r(^ + 1) . . . r(^ + 2=1 



-ds 



q i \ q p> 



q p 



where a = v -l = 2 £ \ n 

9 9 \9. 

We use 



in Lemma [1.21 Then 

/(a,/) 0*0 



2 \ 1 

9 , 



TTm,n 
P,1 



k I (op,i;Ap) 

l(6 9 ,feflg) 



T(gs) x 



pis 



2vrix J L (2 7 r) 1 i £ ^ s -ir(/s)r(/s + ±) • • • F(ls + ^ 



-ds 



(27r)^g 9S -2r(s)r(s + ±) • • • T(s + 2=1) x pZs 



2«x y L (2 7r )^ £ ^-3r(/s)r(/s + i)---r(/ s + 2^1; 



i i^' sr(s)r ( s+ 0'" r ( s+ 



2m 



9-1 
9 



P 



i Zsr(Zs)r(/ s + i)..-r(/ 
i^Tq i / r( s + i)r( s + i)---r( s + 



pi 



x pl q q 



s(2tt)V 2vri 
x{2tt) 2 



9-1 
9 



*r(j a + i)r(fo + ij...r^ + 2=i 

p p/ i(i,o,(|,0.-,(^,0 



ppz 



c/.S' 



a;P'g9l(i,i),(i,i),..,(^i,i) 



Since we have the following 



#9,0 

P,9 



feo= 



a; P yl(i,i),(i,i),-,(^,i)_ 

(-i) fc °r(-fco - g) ■ ■ ■ r(-fcp - f)r(-fco - \] 
\ U r(-ik + ^)T(-ik + a-f ) . . . T(-ik + {p - l l q - plq ) V^V 



p i \ fc +i 



fci=0 

oo 



pq i 



(p-l)q-pl\ \ x P l q1 



pq 



^ + q -^n-h - in-h + $ . . . r(-^ + g 

+ ^ fe! r(-z*2 + E ^)r(-Zfca + *=J*) • • • r(-/A; 2 + V .*-'V 



+ ■ 
+ 

X 

+ 

+ 
+ ■ 
+ 



oo 

V 



-i)*.-r(- Vi + J)r(- Vi - 2 = ? ) • • • r (- Vi - J! 



CT=o Vi! r(-j Vi + ^)r(-^Vi + ^S^) • • • r(-z Vi + 



x pl g q 

x iq 

P 1 



P 



Ip 



x l Pq q J 1 



^ |( -a^,-i)-(-|,-i)(-l,-l) 

-pig i\ f (p-l)g-pig 



X P'g9'(l-/,-0( s 



,-0-(i 

pg 7 ' v pg 



2/2 2 



s*V 1 (»V*,-Q(ffi-0-( (p "^" p ' -0 



(p-l)g-2pi pq-2pl 



x pl q q l(tS r i)-(™,-0( ! 



-0 



we get 



/(a,i)(a:) 



p(«-i)»? 



q-l^p 



(9-2) iw fcj) 



gp' g g | ( pg " ( r i)p ' ,-o( 



g = (g = l)pj 



-!)•••( 



( p _i) g _( 9 _i)pi 



x{2tt) 2 



iW).(i.0.-.( I r 1 .0 



p 



x(2tt 

iyW ^ / p 7 i 



P^ -l)(-i -1) 



x(2tt 



g-g / 



-1* 



x pl q q l(i-i,-o(g^g,-Q-( fe^^g ,-o_ 

p* .(^[^..-Od-iii^-od-iiis'-o-c^-ii],-!'- 1 ) 



1 Xx-'iqi 



a^9« ( a ^-')( s ? ! - , )-( 



(p-l)g-jP' 

P9 ' 



-0 



where 



is the Levy jump function. 



Note that the series are absolutely convergent satisfying the condition X)j=i ~ Sf=i > — 1 

^2 ^2 

-Zp + g - 1 > -1 ^ g > Zp ^ 2 > / ^ 2 > 2 M ^ i >1> (p-Y 1 = a since a < 1. 

r i i r p P q \q P \q J 



If Z belongs to the set {2, 3, 4, • • • }, then we have 



We use 



f(a,i)(x) = 



x(2tt) 



pq _ H i,o 

q-p p,q 



jf l (0,0,(1,0,- ,(^,o 



pq 



x(2ir) 



2ni 



xP i qq \ m ,(i,i),..,(^,i) 



ppi 

x pl q q 



ds 



r(is) = (2vt)^ i ls -^r(s)r i s + i)...r[s + 



i-i 



(2.6) 



Hence we have 

Jpq 1 



x(2tt 



r( a )r( a + |)...r( f 



l (27r)^/ is -ir( s )r (s + i) • • • r ( 8 + 1 -=±) r (is + ± • ■ • r z s 



x pl q q 



fpql 



x(2tt) 2 
y/pql 



r( s + i)...r( s + 



2-1 
g 



x(2tt] 



rrg-1,0 

g+i-p-i - n p+/-2,g-l 



^r( s + f)...r( s + ^)r^ + ij...r(^ + ^ 

^ i(U).-.(¥.i),(^),-.(^^) 



p-l 
P 



Note that for I 6 {2, 3, 4, • • • }, q should be of the form kp, k = 4, 5, 6, • • • and I2 < l\- ^2j=i Bj — 
Y7i=i A-i > means q — 1 — (/ — 1) — l(p — 1) > =>■ q > Ip =>- kp > Ip =>■ k > I. But this condition 
is always satisfied since / = (k) ! i . 



If we put / = 1 in (|2.6|), then 
fpq 



pq 



x{2n) 2 



q-l 



pP (Q,l),(l,l),...,(fci,l) 



rr9 -l,0 



<?-p 
I 2 



E 







y x 3 « g J J 





pp 



x(2tt) 



X(27T, - . =1 

by applying the formula r(/3 + 1 — t> 



' g-jp j\ f (p-l)g-JP _-A 
v pg ' / V pq ' / 



v 1 LgJq-i 



a-iJ-H))^ ^ Lemma 11 



™7 E 



x(2vr 



Q-P 
I 2 



9-1 



>- 
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J 



3=1 



x „-iF ( 



9-2 



■1) 



q-p_ 



*2 = J- \ pg 



i-2=i+m 

<3 LgJq-l 



For the case of / = 1, the condition for their convergence is trivial since p < q. 



■1 



i1 



' p-i 
> p 

'2=1,1 



Some special cases will be given. We will use (2.3) to show some known results. 



Example 2.1 When p = 1, q = 4, a = j, we /iai>e 



/(ii)(z) 
2 



X(27T 
3 



_ 3,0 

3 -'-'0,3 

1 2 



4 % l(i > i),(!,i) 1 (f 1 i) 



E 



x(2n)? \x 3 ±4? 



n 



»1=2 



2xI(2tt)I Vv V2 ' 



-1 

256a; '(!)(!! 



0-^2 



4% l ( 1 -i+[!] 2 )( 1 -!+[l] 3 ) 



8a;f(27r)I V V V " 



256x'(f)( 



6 



1 -r' 1 



1 



32x$ (2n)V V ^ V 2 1 J ' 2 

^G^ r ( 1 + i) r ( 1 + i> 



-l 



-i ._ 

256x'(f)(i 
-1 ._ 
256x '(!)(! 
1 



256a; l (f)( 1 +i; 



256x l ( 1 +l)( 1 +I) 



}• 



We will use (2.2) to show some new results. 



Example 2.2 Let p = 1, q = 2, h = 2, l 2 = 1, a = Z = \/2, Z/ien 



1 ,(l,V2) 



4x v^l(l,D,(|,l) 



V 7 ^ / 1 



4^ 



1*1 



1 (1-V2-V2) 



2x^2 



— 1*1 



1 .(i-i) 
4^ '(l-^.-Vz 



3 A process to generate more representations 



For the same rational a, more representations of one distribution can be generated by using (2.3) 



and (2.4) 



Example 3.1 When p = 1, q = 2, a = | in (2.3), we /lave 



1.0 



-"0,1 



4x l (l- 1 ) 



2^2 1/71 



F 



-1 

Ax 



2^xf XV V^x) 



Example 3.2 Let p = 1, g = 4, £1 = 2, I2 = 1, a = 7 2 = |, I = 2 in (2.4), i/ien we /iave 



/(i,2)(^) 



27TX 



rr3,0 
^1.3 



X 2 4 4l(i,l),(| )1 ),(|,l) 

-1 V 1 



1 ^T(-v+\)(-l 



2nxl £ vl 



4 3 x 2 



4V'i 



4x 



o-^i 



4% 2l (l, 



Note that exp (^) = o-^i 
another representation can 



1 1- 

4 3 x 2 



h o^i 



1 1- 

433,2 I (:p 



and by setting a = (|) 2 = | in (2.4) 
3e born in the sum of faster convergent series. 



4 Levy smashing on the family of gamma density functions 
and Levy-smashed gamma stochastic process 

Mellin transform of a density function in statistics shows its statistical structure and this technique 
can be used as a tool to blend two independently distributed random variables. In this section, 
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we show what Levy effect could be and how we should understand it. To start with, consider 
the Levy density function § L J ~rhzs) x ~ s ds and the gamma density functions Y^e~ x where 



r(i-i) 

v a. a 1 



x~ s ds is the one-sided Levy density function 



< a < 1, < 7 and < x < oo. ^ § L aT (i- s y 
found in j5] constructed in a different way in [2]. We will use the Mellin transformation of the 
form f /ii(|)ft, 2 (t) y, where hx(x) and h%(x) are certain density functions. Then we have 



/(a)( 



;r J 
1 

1 

2~7ri 



1 



r(i- 
«r(i - s] 
r(s + 7 -r 



-x s t s ds 



,dt 



r(7) 



i 

2?ri 



aT(l — s 



-x 



am - s 



r(7) 



x s <is by transformation s = 1 — asi, 



r(si) r(7 -asi 
r(aai) r( 7 ) 

and its Laplace transform 



-x 



<XS\— 1 



oo 

= E 

fc=0 



2ttz 

l) k T(j + ak) 

fc! r( 7 ) 



C?Si 



T(s) r(7 — as) 



OO £S+J-1-1 



~ l dtds 



„as— 1 



r(as) r( 7 ) 



dsdx 



1 

2tu 



(v) 



ak 



r(s)T(7 — as) 



y~ as ds 



(4.1) 



(4.2) 



-y 1 

When a = 1 in (4.1), (4.1) becomes gamma density Y^e~ x . So f(i)(x) is a one-sided function 
concentrated at x — 1 for IR + . To understand its effect, put a = \ in (4.1), then we have 



1 

2vri 



2 si) r( 7 ) 
r(2 S ) r( 7 - s) 



2s 



1 

2t^ 

oo 

E 

fc=0 

4r( 7 + |)(4x)^ 1 



„s-l 



r(«) r( 7 ) 
-i) fc r( 7 + | + 



2ds 



1 

2vri 



fc!v^r( 7 ) 



fc ) -fc-S.-fc-I 
— — X 24 K 2 



2 2s - 1 r( g )r( g + i)r( 7 - g ) 
r(«) r( 7 ) 



r( 7 



;r 



s-l 



2rfs 



2v^Fr( 7 )x 



1 + 



4;r 



1 + 4x)" 



v^r( 7 ) 

Fig. 1 shows the impact on the family of gamma functions. 





gamma family 




Levy smashed gamma family 


(a) 7 = 1 






2(1 + 4x)" 3 / 2 


(b) 7 = 2 


xe~ x 




12z(l + 4x)~ 5 / 2 


(c) 7 = 3 


~X ^ & 


<-> 


60x 2 (l + 4x)- ? / 2 


(d) 7 = 4 


^x 3 e~ x 


<-> 


280x 3 (l + 4x)~ 9 / 2 



(4.3) 
(4.4) 
(4.5) 
(4.6) 
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2* (l+4*t) A ( - (3/2) ) 



exp ( - t ) 



S0*(t A 2)*(l+4*t) A (- (7/2) 




(t A 2*exp(-t) ) /2 



10 12 14 16 18 20 



12 *t* (1+4 *t) A {- (5/2) ) 




Figure 1. 



10 12 14 16 18 20 



(b) 



280*t 3* 
(l+4*t) * (- (9/2) ) 




(t A 3*exp ( -t) ) /& 



6 8 10 12 14 16 18 20 



(d) 



f(a)(x) = 2^ §l Yfcsj F ^(7)^ xas ~ l ds is absolutely convergent series for all x since fi = a — a + l = 
1 > 0, see 0. From the fig. 1, f< a ){ x ) will be called as Levy-smashed gamma density functions 
especially when the parameter < a < 1. Note that a can be any positive real number. 



The stochastic process X(t),t > with X(0) = and having stationary and independent 



increments, where X(t) has the density function ^ 



r(a) T(t-as) „gs-l 



L T(as) T(t) 



X 



l ds, < a < 1, will be 



called Levy-smashed gamma stochastic process. The Levy-smashed gamma stochastic process 
X{t) has the distribution function, for t > 0, < a < 1, F( a)t )(x) = ^ & r(i+<L) F rft) xas ds. 
From pQ and [9], we can prove that the Levy-smashed gamma distribution with parameter a is 
attracted to the stable distribution with exponent a, < a < 1. Namely, 



lim L< 



f*OQ | 

lim / 



jf\ — ) = inn / e 

v n n~s>oo _/ 27u J L T(as) r(n) 



rfs) r(n — as' 



hm — d> v ^n as y 

n-+oo 2m 7 £ 1 (nj 



n as x as 1 dsdx 



ds = ^ytim^(iy k 

n^oo^-^ kl Tin) \nJ 

k=0 K ' 
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5 Remarks 



In [3], they consider the signalling problem for the standard diffusion equation §iu(x, t) = D-^u(x, t) 
with the conditions u(x,0 + ) = x > 0;w(0 + ,t) = h(t), u(+oo, t) = 0, t > 0. And they say 
Then the solution is as follows u(x } t) = J * Qg(x,r)h(t — r)dr, where Qf(x,t) = 2 J^ t~ i exp— -f^ . 
Here Qf(x,t) represents the fundamental solution (or Green function) of the signalling problem, 
since it corresponds to h{t) = 5(t). We note that 

gfa,t)=p LS fcii):=-j^exp(-^), t>0, /i=^ (5.1) 



where PLs(t; p) denotes the one-sided Levy-Smirnov pdf spread out over all non-negative t (the 
time variable). The Levy-Smirnov pdf has all moments of integer order infinite, since it decays 

3 

at infinity as t~%. However, we note that the absolute moments of real order v are finite only if 
< v < |. In particular, for this pdf the mean is infinite, for which we can take the median as 
expected value. From VLs(t me d', aO = §, it turns out that t me d ~ 2/x, since the complementary 
error function gets the value \ as its argument is approximatively 



l 

2 ' 



With the inspiration from the above paragraph, take the Levy density function with the pa- 
rameter a — then the density function is well known to be ^_ 3 exp{— ^). Put this in the mellin 

convolution formula f. /ii(|)/i 2 (£)y, then it becomes J °° 3 exp (— 4-) hzifydt. 3 exp {— jA 



has the same form with (5.1 ) when t — p. Therefore the cases of Levy smashing treated in section 
4 can be thought of as superstatistics in Physics and Bayesian statistical analysis, subordination 
in statistics, namely, 

Jo 2^x1 t Q k\ (4x)* T(t) e dy 2 ^V(t)xi t Q k\ (4x)" J V 
1 (~l) fc r(t + k + \) _ 4T(t + \)(Ax)^ x 

2y/*T{t)x% kl W " ^(t) 1 + X) 

But in this paper, the concept of Levy smashing is totally different from the point of view of 
superstatistics in Physics and Bayesian statistical analysis, subordination in statistics. 
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